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We study the T-violating lepton transverse polarization (Pf 1 ) in three 
body semileptonic heavy meson decays to pseudoscalar mesons and to vector 
mesons. We calculate these polarizations in the heavy quark effective limit, 
which simplifies the expressions considerably. After examining constraints 
from CP conserving (including b — > 57) and CP violating processes, we 
find that in B decays, P of the muon in multi-Higgs doublet models can 
be of order 10%, while P 1 - of the r can even approach unity. In contrast, 
P^ in D decays is at most 1.5%. We discuss possibilities for detection of 
Pi~ at current and future B factories. We also show that in decays to 
vector mesons, unlike in decays to pseudoscalars, can get contributions from 
left-right models. Unfortunately, Pf~ in that case is proportional to Wl-Wr 
mixing, and is thus small. 
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1 Introduction 



The Standard Model (SM) has thus far met with incredible experimental 
success. Nevertheless, many hypothetical extensions to the SM remain phe- 
nomenologically viable. Since new physics often provides new sources of CP 
violation (CPV), one good way to search for such extensions is to consider 
CP violating observables which are negligible in the SM, but which can have 
large contributions from other sources of CPV. 

A major barrier to any candidate for such an observable is the upper 
bound on the electric dipole moment of the neutron, d n , which is now 
around 10~ 25 e cm The SM explanation for CPV, the Cabbibo-Kobyashi- 
Maskawa (CKM) mechanism has come to be accepted by many as the 
source of CP violation in the neutral K sector not only because it predicts e 
to be in the right range, but also because it predicts d n to be negligible M. 
As the upper bound on d n has plummeted, many potential explanations for 
e from other sources have run aground, and thus it is more difficult to find 
observables which have good prospects of detecting CPV beyond the SM. 

One such observable is the transverse polarization of the lepton in semilep- 
tonic K^ 3 decays [|J, P^, which is the T- violating projection of the lepton 
spin onto the normal of the decay plane, i.e. P z ~ si ■ (k x p) ||, where k 
and p are decay product momenta. It arises from the interference between 
two amplitudes with non-zero relative phase. In practice, one measures the 
asymmetry between the number of particles parallel and anti-parallel to the 
normal of the decay plane, 
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pt = — — -. (i) 

There are several advantages to using such a semileptonic CP violating ob- 
servable. First, semileptonic decays occur through a single SM diagram at 
tree level, so P t is negligible in the SM pi. Thus a non-zero P t is a sig- 
nal for new physics. Second, theoretical uncertainties in semileptonic decays 
are much smaller than in purely hadronic decays. Finally, P^ in semilep- 
tonic decays comes from both the quark and lepton sectors, so that purely 
hadronic or purely leptonic CP violating observables, such as d n or d e , do not 
necessarily strongly constrain P^ [[F]. In fact, there exist reasonable models 
for which P t in decays can be of order 10~ 2 -10~ 3 , consistent with all 
other constraints || [|. Such values are well within reach of experiments. 
The last measurements of P t were done at Brookhaven National Lab on 
K + — > 7r°// +l/ /i decays. Their combined result was [Kj 

P^ (K + -> 7r°/i + ^) = -1-85 ± 3.60 x KT 3 , (2) 
which implies a 95% confidence upper bound of about 0.7%. There is also an 



experiment currently under construction at KEK |LT] which hopes to push 
this bound down by a factor of ten ||12|| . 

In this paper we consider P t in heavy meson decays of the type M — > 
m I vi and M — > m* I z//, where M and m are pseudoscalar mesons, and m* is 



a vector meson. P t has been studied in decays to pseudoscalars JT3|, but 
not in decays to vector mesons. We derive expressions for in M — > m* I v% 
decays, as well as in M — > m I V\ decays, in the heavy quark effective limit. 
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This greatly simplifies our results. One can even obtain analytic expressions 
for Pj 1 , the polarization averaged over all kinematical variables. 

In decays to pseudoscalars (M — > mli/i), P t is sensitive only to spin 
effective Lagrangians j|, |HJ, which makes it a good tool for probing non-SM 
Higgs physics ||. We find that this holds for P t in decays to longitudinally 
polarized vector mesons (M — > m* L I i>i), but that P^~ in decays to transversely 
polarized vector mesons (M — > / v\ and M — > m^ 2 I v\) is sensitive only to 
new V and A physics, such as left-right models. Unfortunately, Pf- in that 
case is proportional to Wl-Wr mixing, which is constrained to be small. 
However, in the former case, multi-Higgs doublet models yield encouraging 
results, even after imposing CP conserving and CP violating constraints. 
There are reasonable models in which P^~ in B — > Dtv decays can even 
approach unity. 

Section 2 lists the form factors needed for our calculation. We consider 
contributions to P^~ from multi-Higgs doublet models in Section 3, and from 
left-right models in Section 4. Possibilities for detecting P^ in various decay 
modes are discussed in Section 5. 

2 Form Factors 

From Lorentz invariance and basic symmetry considerations, we can write 
the hadron matrix elements (HME) for decays of pseudoscalar mesons (M) 
to pseudoscalar (m) and vector mesons (m*) as 
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(m(fc)| V, \M{K)) = f + (K + k), + f-{K- k), , 



(m(k)\A fl \M(K)) = 0, 

~M 



(m*(k,e* x )\ V, \M{K)) = Li (e^e^K^ 



A 



(m*(k,e* x )\ A, \M(K)) = A x Me* x + -f • K)(K + k) 



M 



^3 



+-^(ej.lf)(if-fc)^ (3) 

where for M + decay, = Dj^U and A^ = D r y l /y 5 U (U and D are the 
appropriate up- and down- type quarks for M and m). The axial vector 
HME for M — > m is zero because there is no way to form an axial vector 
with just K a and k^ . We have used M and m to represent both a meson and 
its mass. The form factors f±, V x , and A x _ 3 are functions of (K ■ k) and r = 
m/M. Here A is the polarization index. We will refer to the m* longitudinal 
polarization by the label A = L, and the two transverse polarizations by the 
label A = Ti, T 2 , for e\ in the decay plane and perpendicular to the decay 
plane, respectively. 

From these vector and axial vector HME, one can derive scalar and pseu- 



doscalar HME using the Dirac equation [13 



(m(k)\ S \M{K)) = - r 2 ) + f_t) , 

(m(k)\P\M(K)} = 0, 

(m*(k,e* x )\S\M(K)) = 0, 

(m*(k,e* x )\P \M{K)) = ~ M (e* x ■ K){A X + A 2 (l - r 2 ) + A 3 t),(4) 
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where for M+ decays, S = D U, and P = D'fU, and t = (K - kf/M 2 . The 
masses {mr,,mu) are (mj, m c ) in B decays and (m s ,m c ) in D decays. The 
middle two parity-odd matrix elements in (f|) are zero because there is no 
way to form a pseudoscalar using only K a , k 13 and ex 1 . Note that the factor 
{e\ ■ K) implies that (m*(k, e* x ) \ P \M(K)) is non-zero only for longitudinally 
polarized vector mesons. 

Recently there has been a lot of interest in heavy quark effective theory 
(HQET), which considers the limit M, m — > oo. The principle tenet of 
HQET is that (v'), the four- velocity of M (m^*'), is unchanged by QCD 
corrections [IB]. Thus it makes sense to write the HMEs in terms of velocity 



(m(y')\V tt \M(v)) = VMrn + v'Y + - t/) m ) , 
(m(v')|A M |M(v)> = 0, 

^ |M(i;)) = zVM^£ Vl (e^eJVV) , 
(m*« s* x )\A, \M(v)) = VM^(ua + v ■ t/)e^ - Ufo ■ «K 

-U(eA'«H), (5) 



From (R|) and (pf), one can derive relations between the form factors 18 



k,= — U„, -»-4=«- 



^i = — — ► — t^S, 



6 



* = +^(U-rU)-+^, (6) 

where x = (K ■ k)/M 2 = rv ■ v'. In the M, m — > oo limit, £ + = ^ = ^ = 
£a 3 = £ and £_ = £a 2 = 0) so that all the form factors can be written in 
terms of the Isgur and Wise function, ||19|| . Note that the HMEs are 
normalized so is equal to 1 at zero recoil (x = r or v ■ v' = 1) [pOj] . 
Specific forms for are listed in the Appendix. 



3 Higgs models 

3.1 Transverse Polarization 

As we said, semileptonic pseudoscalar decays to pseudoscalar mesons, M — > 
m I u, and to longitudinally polarized vector mesons, M — > m* L I u, can arise 
only from the interference of new scalar physics with the SM. In this section, 
we consider contributions to P L (M — ► mlvi) and P ± (M — > m* L lui) from 
models with charged Higgs scalars. Other types of contributions are possible, 
such as from scalar leptoquarks J|, §. 

T. D. Lee first proposed that CP could be violated via phases in a model 
with two Higgs doublets |2l| . This idea was refined by S. Weinberg with 
the elimination of flavor changing neutral currents (FCNC) in a model with 
three Higgs doublets, using a symmetry to ensure that only one Higgs doublet 
couples to each right-handed fermion field; what is commonly referred to as 
natural flavor conservation (NFC) |22j . There are various other ways to avoid 
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the FCNC problem p3| , |24|| , but for simplicity, we concentrate on models 
where NFC is either exact, or partially broken f25fl . We will assume that the 
CKM phase is non-zero, so we do not impose strong constraints on CPV in 
the Higgs sector from e. Even if CP is broken only spontaneously, a non-zero 
CKM phase can arise after integrating out super-heavy fields, so we see no 
reason to take it zero. 

We are interested in the interference of a charged Higgs boson with the 
SM W boson, so one need only parametrize an effective Lagrangian for the 
charged Higgs coupling to fermions. In a model with N charged scalar fields, 
one obtains a Lagrangian in terms of the N — l physical charged Higgs bosons 

& 

Y-l. 

V 



i JV-1 

-C H+ = - £ \u t U L V L M D D R H+ + fiJJjiMuViD^Hf 



+ li N L M E E R Ht]+H.c. (7) 

Here v is the SM Higgs VEV, v = (4G F / V / 2)" 1/2 ~ 174GeV; U, D, N, and 
E are fields for the up quarks (U T = (u c £)), down quarks, neutrinos and 
charged leptons; Mo, Mjj and Me are the diagonal mass matrices; and Vl 
is the CKM matrix. 

If the coefficients ati, and 7« are complex, the interference between 
the charged Higgs and W boson amplitudes in Fig. 1 produces a T- violating 
transverse polarization of the lepton. Since the H + amplitude is proportional 
to the matrix elements in (^), one gets contributions to P l (M m*lv) only 
for decays in which the m* is longitudinally polarized. This means that if one 
can veto decays with transversely polarized m*'s, the denominator in (HD will 
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be reduced while the numerator will remain unchanged, leading to a larger 
polarization. 

Let us evaluate P^ in terms of the velocity dependent form factors. Then 
we can take the heavy quark effective limit, which allows us to write with 
only one form factor, £(x). We calculate P t for semileptonic pseudoscalar 
decays to pseudoscalar mesons, to longitudinally polarized vector mesons, 
and to unpolarized vector mesons in this limit: 



H+ 



3vr (1 - r 2 )(x + r)(x 2 - r 2 )^t £(x) 



2 



™ l*--*«*u)= C H+ - ± ^, (8) 



h+ * r \ ^ 3tt (1 - r 2 )(x + r)(x 2 - r 2 )y/i f( 



x 



,2 



P^x) ( M ^ m* L lu) = C H+ — ± n "J,,.", 7V ttzt 2 , (9) 



4 (1 — r) 2 (x + r) 2 Xi £(x) 



H+ 



P^lx) M ^ m*lu ) =C H+ — 



3tt (1 - r 2 )(x + r)(x 2 - r 2 )^/i £(xf 



4 (1 — r) 2 (x + r) 2 xi + 4t(x + r)xx\ £(x) 2 

(10) 

We list the full expressions with general form factors in the Appendix. Note 
that we have already integrated Pf 1 over one kinematical variable ((K-p)/M 2 ) 
so that Pf 1 is only a function of the remaining kinematical variable x (where 
x = (K ■ k)/M 2 and X\ = \J x 2 — r 2 ). This integration gives the factor 3tt/A 
in For M+ -> m 0{ *H + ui and M° -> m~^l + Ui decays, the new 

physics coefficient is given by: 



Mm/ * 1 / m D m v \ 

^+ = 172- L T72- ( — — — lma ^ + — — -~ Im A7i , (11) 
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while Ch+ for the CP conjugate decays has the opposite sign [P^j . Here 
mi, mjj and rri£, are the lepton and current quark masses specific to each 
decay, and Mr. and the coefficients aa, fa and 7, come from the effective 
Lagrangian (|7|). The — (+) applies to M — > mlv [M — > m*lv) decays. Since 
mu > in D decays, it follows that P l ± (D + —>■ K°l + i/) has the opposite 
sign as P[ L (K+ -> tt°1 + u), P l ± (B+ -> D°^l + u), and P Z ± ( J D + -> K°*l + v). 
It also means that Cy+ in the decays to m* are somewhat suppressed over 
those to m when and my are of the same order, as in B decays. 

We have neglected all lepton mass effects in the denominator of (3)- fll0|) ■ 
For I = fi, this is always a very small effect. In I = r decays, it changes our 
results only qualitatively when P 7 1 ~ 1, i.e., when H + effects are important in 
the denominator of flU). In that case, it might be possible to see new physics 
effects in changes to the branching ratio of B — > D^*'tu. 

To estimate the size of P^ in various models, we must integrate over the 
remaining kinematical variable x. In an experiment, one generally measures 
the overall asymmetry in ([I]), rather than measuring Pi (x) for each x and 
then averaging. So we must integrate the numerator and denominator of 



)-(10) separately: 



J? (M 5 mil.) = C„+f£, 
iF(* 5 „•!„)= C^jJ^-, (12) 
where I±, Is, II an d It are integrals of the kinematics in (^)-(^). Unfor- 
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tunately, this means we must know something about the overall form factor 
£(x). In the Appendix, we list two possible parametrizations for £(x): a 
relativistic oscillator model, and a monopole approximation. P^- 1 in decays 
to pseudoscalars in these models differs by at most 15% for r in the region 
of interest (r > 0.25 for all the decays we study), and considerably less for 
decays to vector mesons. If we set the monopole parameter p equal to 1 in 
(pip, we can obtain analytic expressions for in terms of r. We list the 
corresponding Ps in the Appendix. From Fig. 2, we see that choosing p = 1 
instead of 1.2 (in order to obtain analytic expressions) changes P^ by only a 
few percent (for r > 0.25). Even naively dividing out £,(x) 2 from the numer- 
ator and denominator of (^)-([T0|) gives results which (for r > 0.25) differ by 
30%, or much less, from the other parametrizations of 



3.2 Constraints 

For the purposes of placing constraints on Pf 1 , we make two simplifying as- 
sumptions. First, we take the «j, $ and 7$ to be flavor diagonal. This 
strictly holds only in models with NFC, so Higgs models without NFC may 



have somewhat weaker, more model-dependent bounds [27]. Second, we will 
assume that the lightest charged Higgs mass eigenstate, h + , gives the domi- 
nant contribution, so that we can drop the subscript i on the coefficients a, 
(5 and 7. In 3HDMs, Imo;i7* = — Im«272 an d I m /3i7! = — I m /^272 5 so m that 
case we are simply making the replacement M~+ — M~+ — > M^+. This has 
virtually no effect on CP violating constraints, because they have the same 
behavior, and the CP conserving constraints tend to require a large splitting 
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between M H + and M H + anyway. 

We want to constrain C H +, which now depends upon Imary*, Im/^7*, 
M^/Ml+, and the masses involved with M and rrS-*\ In the general case 
(given our two assumptions), we can bound Imaj* M^/M% + directly from 
the experimental upper bound on P^(K + — > 7r°yU + z/ M ) of 0.7% to obtain 

A/f 2 

|Ima 7 *| -f < 730. (13) 



Since rrtu is small, P^(K + — > ir /i + z/ M ) is insensitive to Inu37*. The best we 
can do is to use |Im/37*| < • |7|. From the bounds placed upon and | 7 | 
by [p8 |, we obtain 



|Im/? 7 *| < 160^ < 285, (14) 

1 1 M 2 + M /i+ v y 

From flnp, one sees that the upper bound on Im/37* My V /M^ + decreases with 
increasing M h + and is at its maximum when M h + is at the model independent 
lower bound of Mz/2. We can use ( |13"D and ([[4]) in (|l"2"D to obtain upper 
bounds on P^ for various decays. Our results are summarized in the first 
column of Table 1. 

Let us now specialize to the case of 3HDMs. The CP violating coefficients 
can be written IMI 



Im«7* = i sin 26 3 sin 5^-, 
Im/37* = ± sin 2fl 3 sin <5^, 

Ima/F = | sin 2fl 3 sin 6-^-, (15) 
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where v u , va, and v e give mass to the up quarks, down quarks and charged 
leptons, respectively. 9 3 (5) is a free, CP conserving (CP violating) param- 
eter of the model. For convenience, let us define 

M 2 

k = Isin2# 3 sin 51 — ^, (16) 

I I M 2 + > V ! 

so that Iincry* My V /Mf L+ and Im/^7* M^/M|+ are just given in terms of k 
and the VEVs. The relations in ( |i5|) are not enough by themselves to better 
the constraints given by (|0|) and flR]) , so we consider specific models. 

A common assumption is that the three VEVs are all of the same order, 
i.e., v u ~ Vd — v e . We refer to this as the VEV Equality (VE) model. In the 
VE model, all three CP violating coefficients are of order one, and P^ will be 
quite small. But with one VEV for each type of massive fermion, this need 
not be the case. Since fermion masses are proportional to the VEVs as well 
as the Yukawa couplings, it is quite reasonable to suppose that the hierarchy 
in the fermion masses lies in the VEVs, and not the Yukawa couplings M. 
Suppose the third family Yukawa couplings are of the same order. Then one 
has v u : Vd '■ v e ~ m t : m& : m T , which implies that 

lwl w~ ~ ^r K - (17) 

M h+ rribm T 

so that Pf 1 need not be small ||. We will refer to this as the VEV Hierarchy 
(VH) model. While the VH model provides a reasonable justification for 
considering large ratios of VEVs, it does not solve all the mass hierarchy 
problems. We view the VE and VH models as two reasonable extremes, 
much in the same way that the range 1 to m t /mb is considered for "tan/3" 
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in 2HDMs. 

For simplicity, we define the VEVs in the VE model to be identically 
equal, and in the VH model to have the ratio m t : : m T exactly Since 
fllDD implies k < 3.2 or so, P^(K + — > 7r /i + z/ M ) does not put any further 
constraints on k in the VE model. However, the VH model can reach the 
upper bound on P^(K + — > 7r°/i +l/ £t)> and one needs k < 0.5. We now must 
consider if there are any other constraints on k which would force P/~ to be 
small. 

As we said in the Introduction, the most stringent constraint on CPV 
often comes from the electric dipole moment of the neutron, d n . The purely 
hadronic coefficient Ima/3* is very constrained by d n |§ , and in the VE model 
we find that we need k < 1.2. However, in the VH model, Imcry* /Ima/3* is 
large, and the upper bound on k is only about 6, which is ten times weaker 
than the bound. This is a consequence of the semileptonic decay — only 
quark-lepton CPV is enhanced in the VH model. 

CP conserving processes may also constrain P^. Consider the inclusive 
decay b — > 57, whose branching ratio is now bounded below 5.4 x 10~ 4 at 



the 95% confidence level by the CLEO collaboration [29|. This FCNC decay 
occurs via a one loop diagram in the SM with a branching ratio of 3-4 x 10 -4 
[ j30|j . In 2HDMs, the charged Higgs contribution adds constructively with the 
SM contribution, and one can bound the charged Higgs mass to be above 



about 450 GeV []3T | . One would like to generalize this result to 3HDMs. The 
amplitude for b — > 57 (at the W mass scale) can be written |32 
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Nik) ' (18) 

where the sum over % runs from 1 to N—l, and one can show that /3j/3* = (v 2 — 
w 2 )/f 2 , Re«i/3* = 1, and Imaci/3* = 0. In the SM, only the first term is non- 
zero. For iV = 2, we recover the 2HDM limit, i.e. (|/?| 2 , Rea/?*, Ima/3*) -> 
(v^/v^,l,0). In 3HDMs, one can have cancellations between the pieces as 
long as Hi and H% are not degenerate in mass. It turns out that for both 
the VE and VH models, Reaifil can be less than zero, so that for sufficiently 
large M H +, there is no bound from B(b — > s'j) on M H +. For M h + ~ My/ 



(or smaller), sin 29% sin 8 must be somewhat smaller than one [33], but this 
is not enough to better the constraints on k we have derived thus far. 

There are also constraints from B(b — » 57) on Ima/5* p2| , [35[| , which in 
turn constrains Pj 1 via (^). Since the last term in ( |IBD is purely imaginary, it 
does not destructively interfere with the other terms, so that the contribution 
from lmaj3* to B(b —>■ 57) is always positive. However, even for ~ 



Mz/2, one can only bound Ima/3* < 2 ||35|| , which is satisfied in both the 
VE and VH models. Since the CLEO observation of B(B -> K*-f) = (4.5 ± 
1.5 ± 0.9) x 10 -5 |2^| effectively sets a lower limit on B(b — > 57) of about 



10~ 4 , the constraint on Ima/3* from b — > S7 will never be able to strongly 
constrain P t in these models. 

Finally, we note that the VE and VH models give specific predictions 



for Im/^7* (see (13)), and in both cases it must be less than 2. In general 
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3HDMs, one cannot improve upon the bound in ([14]), though large Im/37* 
would require small v u /vd as well as very large v/v e , which is not as appealing 
theoretically as either the VE or VH models. 

In Table 1, we summarize the maximum values for P^ allowed in the 
VE (VH) model, with a bound of k < 1.2 (/c < 0.5) coming from the upper 
bound on d n (P^ in K decays). 



4 Left-Right Models 

Decays to vector mesons, M — > m*lv, have one more 4-vector than M 
mlv decays with which to construct hadronic matrix elements. The m* 
polarization vector lets us construct both a vector and an axial vector current 
(see (g|)), allowing a non-zero V and A interference term. The upshot is 
that Pi (M — > m*lv) gets contributions from spin 1 effective CP violating 
Lagrangians as well as those of spin 0. 

Let us therefore consider left-right models PEfl , whose charged gauge bo- 
son couplings to fermions can be parametrized by the following effective 
Lagrangian: 



~ Cw+ ~ T2 



9R 
V2 



U Rlll V R D R \W^ + H.c. } (19) 

where Vr is the right-handed CKM matrix. We neglect right-handed currents 
coupled to leptons because they yield polarizations proportional to m„. This 
means that P^ must arise from the interference of the SM Wl diagram and 
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a diagram containing Wl-Wr mixing (see Fig. 1) [|J. We define the mixing 
angle Q by 



) 









(20) 



t W 2 J V - sin C cos C, 
where Wi and W2 are the two mass eigenstates. The interference between V 

and A HMEs vanishes for longitudinally polarized m*'s, so Pj 1 is only non- 
zero for transversely polarized m*'s. Further, the numerator of ([!]) has the 
same magnitude, but opposite sign, for m*'s with Tl and T2 polarizations. 
Therefore, the polarization in the sum of decays to both transversely polar- 
ized m* states, P; _L (M — > m* T lv\ is identically zero, and we must consider 
Pi~ for either Tl or T2. We again write Pj 1 in the heavy quark effective 
limit, 



w+ \ _ 3tt (x + r){x 2 - r 2 )y/i £(x) 2 



P l L (x)[M -4 m* TX lv\ = a 



W R 4 2t{x + r)x 1 {x-r/2) £(x) 2 



PHx)lM W -$m* lv\ - -C + - ( x + r )( x2 - r2 )^ Z(*) 2 (21) 
P x \x)\M — > m T2 /z/j - G w + 4 2t(x + r)xi(x + r/2) i[x? W 



and list the full expressions in the Appendix. The coefficient, 

C w ;^Untl m (^f), ,22) 

depends upon the W^-Wr mixing angle (, the left and right CKM elements 
Vl 3 and Vr = U,D), and gauge coupling constants and qr. 

We can find an averaged polarization by integrating the numerator and 



denominator of ([H]) over x 
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R 4 I T1 

px(m3^)= _ C< |Mizl!). (23 ) 

We again use the p = 1 monopole expression for which results in the 
Iti and It2 listed in the Appendix. We have normalized the J's so that 
In + It2 — It- Fig- 3 shows that using p = 1 (to obtain an analytic 
expression) instead of 1.2 is a good approximation since £(x) 2 appears in 
both the numerator and denominator in (|2"T|). 

Let us consider constraints on P^~ in B T 3 decays. Our Lagrangian in (|7|) 
gives a tree /ewe/ contribution to e' |37]], and we can relate and e'. If 
Im(V^ rD /V L c/D ) is roughly the same order for all UD, then P~F ~ 10~ 2 (e7e), 
which is tiny. It is in principle possible that lm{V^ d /V^ d ) ~ Im(V^ s /V^ s ) ~ 
while Im(V^ /V£ ) ~ 1, which gives P^~ ~ 2£. Nevertheless, \(\ is con- 
strained to be less than about 6% from decays |38 |, and less than about 
2% from b — > 57 |39|], so that we can bound P^- to be less than about 4%. 



5 Discussion 

Let us consider the various decay modes. In particular, we discuss whether 
one should study charged or neutral decays, of B or D mesons, to pseu- 
doscalar or vector mesons, with I — p or I — r. 

Technically, the transverse polarization, P^, is motion reversal violating, 
which is equivalent to T violation only in the absence of final state interaction 
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(FSI) effects [EO]. This is irrelevant in charged decays, e.g. M + — > m°l + ui, 
because they have only one charged decay product, and FSIs are negligible. 
In neutral decays, e.g. M° — > m~l + i>i, there are two charged particles in the 
final state, so one can expect FSI effects of order olemI^ PI| - For this reason, 
measurements of in decays are done on the K + — > -k q ^l + v^ mode. 
But if the experimental sensitivity to P/~ in a given decay is only at the 
percent level, one can study decays of neutral mesons as well. Actually, since 
both B and D mesons are produced in pairs, one must be able to determine 
the charge of the lepton (because P^ flips sign for the CP conjugate decay) 
so that one effectively measures the asymmetry 



1 r 



Acpv = 2 



Pt(M ^ml + ui) - F,(M -> ml'Pi) , (24) 



which is a true CP violating observable. Since FSI effects cancel in Aqpv, 
charged decays are in principle not preferable to neutral decays. 

From Table 1, it is clear that B decays give larger P^ than D decays. 
One can see from (0) that this has two causes: Mo is smaller than M#, and 
the heavier quark mass in D decays, m c , is proportional to Im/?7* instead 
of Im«7*. The former coefficient is more constrained than the latter, and 
models in which Im/?7* is large tend to be less theoretically appealing. For 
example, in 3HDMs, one would need v u /vd to be small while v/v e is very 
large. 

Let us estimate the number of decays necessary to see a 5a signal of Pf 1 
with the maximum allowed values in the general case (column 1 of Table 
1). We use N = 25k/P^ 2 , and take k ~ 10. One needs about 2.5 x 10 4 
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B — > Dfiu decays, which (naively) translates into 2 x 10 6 B's (including B ± , 
B°, and B°). In B — > D*\iv decays, P^~ is 6 times smaller (if one does not 
veto decays to transversely polarized m*'s), but the branching ratio is about 
3 times larger, so one needs about 12 times as many 5's as in B —>■ Dfiu 
decays to see a signal. By contrast, one needs about 1.3 x 10 6 (2.8 x 10 7 ) 
D — > Kfiu decays to observe P^~ for the maximum value in the general case 
(VH model), which naively requires 4 x 10 7 (10 9 ) .D's. 

To observe P^~ of a muon, one needs to stop the muon so it can decay. At a 
symmetric B factory, such as CESR or DORIS II, the muon in B — > Dfiv will 
have momentum of up to 2.3 GeV, which would require perhaps 1.3kg /cm 2 
of material (e.g. ~ 4.5m of Al) to stop it [42]. Unfortunately, even 2.5 x 10 4 



B — > Duu decays is probably out of reach of either machine. Stopping muons 
would be more difficult at the asymmetric SLAC B factory, since the muon 
momenta will be higher in the lab frame, but if it could be accomplished, the 
luminosity should be sufficient to see a 10% polarization. One could consider 
measuring P^ at a hadron collider, where the number of B^ and D^ 3 decays 
would be much greater, but the hurdle of stopping the muon would need to 
be overcome. 

A better possibility may be B t3 decays, because one can have P^- ~ 1. 
One needs perhaps 250 B — > Dtv decays, and 3000 B — > D*tv decays to see 
a 5<t signal. Both of these are probably out of reach of the existing symmetric 
machines, but should be no problem for the SLAC B factory. Unlike muons, 
taus do not need to be stopped, and one can measure the polarization of the 
r from its decay spectrum |j43"| . In r — > n^Ur decays, for example, the decay 
width has the behavior dY ~ 1=(=P T ± ~ 1 — P^r+ cos 9 |4"4]| , where P T ± is the 
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polarization vector of the r , p w is a unit vector in the pion direction, and 9 
is the angle of p n from the normal of the B decay plane. The main problem 
with B t3 decays at the SLAC B factory may lie in defining the decay plane, 
since the 5's do not decay at rest, in which case we may have underestimated 



k 10, 11 



Finally, we note that P^ from left-right models is probably unobservable 
at the SLAC B factory. In addition to the small values for P t required by 
the bounds on Wl-Wr mixing, one needs to measure the polarization of m* 
as well as of I, so that our k is perhaps 100 or more. For P^- ~ 4%, one needs 
more than 10 6 B — > D*rv decays. 

We have derived expressions for the transverse polarization of the lepton 
in semileptonic meson decays, in the heavy quark effective limit. Reasonable 
multi-Higgs models can give a muon polarization in B decays of order 10% 
and a tau polarization of order unity. Both of these should be within the 
luminosity reach of the SLAC B factory, though the tau polarization has the 
advantage of not requiring a stopper. Should a non-zero signal be observed, 
implying the existence of physics beyond the Standard Model, the best place 
to study P^~ would be at a high luminosity symmetric B factory. 

Acknowledgments 

I sincerely appreciate helpful conversations with G. Kane, J. Bernabeu, J. 
N. Ng, D. Bryman and Y. Kuno. This work was supported in part by a grant 
from the Natural Science and Engineering Research Council of Canada. 



21 



Table 1 



Decay 


general 

case 


Pi 1 
VE 

model 


Pi- 
VH 

model 


K + -> 7T° u + v a 


0.7% 


1 x 10~ 5 


0.7% 


D + -> i?° u+ za. 


1.4% 


4 x 10~ 5 


0.3% 




0.51% 


2 x 10~ 5 


0.1% 


D+ -> //+ ^ 


0.27% 


1 x 10~ 5 


0.06% 


B + _> Z)° M + ^ 


10% 


2 x 10" 4 


9.7% 




3.3% 


6 x 10~ 5 


2.9% 


B+ -> D°* ^+ ^ 


1.7% 


3 x 10~ 5 


1.5% 


5+ -> D° r+ i/ T 


~ 1 


3 x lO" 3 


~ 1 


5+ -> £)°* r+ i/ T 


~ 55% 


1 x 10-3 


~ 50% 


5+ -> D°* r+ i/ T 


~ 29% 


0.5 x 10~ 3 


~25% 



Table 1: Maximum values of the transverse polarization {P l ) for 
various decay modes due to SM interference with charged Higgs 
bosons in the general case and in two specific models. K®* {D°l) 
refers to longitudinally polarized K°*'s (D°*'s). For the VE (VH) 
model, we use k < 1.2 (k < 0.5) as derived in the text. Numbers 
of order one are approximate since we neglect H + effects in the 
denominator of ([!]). 
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Appendix 

So far, we have used the heavy quark effective limit, in which £+(x) = 

£vi(x) — £ai( x ) — £a 3 ( x ) — £,( x ), an d €-{ x ) — £,A 2 ( X ) = 0- F° r completeness, 
we list the expressions for the polarization without that simplification: 

V / 4 x'l 

[(l + r)e + (x) + (l-r)C-(x)] x 

[(l+r)Z + (x) + (l-r)U x )} 2 
[(1 - r )(x + r)£ + ( x ) + (1 + r )( x - r)£_(x)] , (25) 



V / 4 (x + rj 2 xi 

[UA X )( X + r)(x- r 2 ) - (U 3 ( x ) + r£ A2 {x))x\] 
- r 2 ) - (^(z) + r£ A2 (x))(x - r)J 

£ Al (x)(z + r) - -(U(^) +^a 2 (z))(1 -r 2 ) + 



(26) 



Pt{x) (M H -X m'lu) = C H+ ^ 
' V / 4 (x + r)xi 

^ Al (x)(x + r)(x - r 2 ) - (^(ar) + rU 2 ( x )) x 

2 



X 



a x (z)(x + r) - + rUW)(l - r2 ) + - <A 2 (z))t 

r 2 i 2 

(x + r) ^ Al (x)(x - r ) - (£4 3 (x) + r£4 2 (x))(x - r)J + 

2r 2 t [^(x) 2 (x - r) + £ Al (x) 2 (x + r)])"' (27) 

where is given in (|TT|) . The corresponding expressions for LR contribu- 



tions are: 
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Pl \ X ) ( M W J m* T1 lv) = ^ + ^ (x + r) ^ : 
^Ai(»)^Vi(ar) 



|£v\(x) 2 (:r - r) + |^ Al (a:) 2 (x + r) 



P^x) M -4 m * T2 lu 



37r (x + r)x\^/t 
w « 4 2t(x + r)xi 



(28) 



(29) 



J^(x) 2 (a; - r) + |^(x) 2 (x + r) 
where C w + is given by ([22]). 

To find the average polarization, we must integrate both numerator and 
denominator over x. For this, one must know £(x). One possible choice 



comes from a relativistic oscillator model 18 



2r 



(x + r) 



g (x—r) 



(30) 



where (3 ~ 1.85 ||18[| . Another possibility is a monopole approximation, 



(31) 



1 + p 2 (x — r)/r 

where p ~ 1.2 ± 0.25 f45f . For most choices of the integration over a; 
must be done numerically, but for the monopole approximation with p — 1, 
one can obtain reasonably simple analytic expressions (see below). Since 
£,(x) 2 appears both in the numerator and denominator of (H), P/~ is fairly 
insensitive to the choice of £(x). We find that for decays with the lowest 
value of r (~ 0.25), the difference between Pj using £(x) from (|3~ID for p = 1 
(analytic case) and p = 1.2, and (|30|) is no more than 15%, and considerably 



24 



less in most cases (see Fig. 2 and 3). Thus we use (|5ID with p = 1 to obtain 
the following analytic expressions for the integrals in ([12]) and (|23|) : 



(l+r 2 )/2 



dx (1 - r 2 )(x + r)xlVt^(xf 



— r 2 (l — r)(l + 6r — 6r 5 — r 6 
15 V A 



2r(l — r + r 3 — r 4 " 1 

7TT 



tan 1 



1 — r 
Vl + r 



(32) 



(l+r 2 )/2 



da; (1 + r)X£(x) 2 

lr 2 (l-r)(l + r) 3 (l + 10r 2 + r 4 ) 3 4/l , 2l 

+ -r 4 (l + rr lnr 

8 1 + r 2 2 K J 



(33) 



(l+r 2 )/2 



c?x (1 — r) 2 (x + r) 2 xi £(x) 2 
1 r 2 (l - r) 3 (l + r)(l + 8r - 6r 2 + 8r 3 + r 4 ) 



2r 4 (l -r) 



tan 1 



1 + r 2 
2r 



1 - r 2 



71 1 

in r 

2 4 



(34) 



'Tl 



Kl+r 2 )/2 

It = j dx 4t(x + r)xxi^(x) 

Jr 



6 



r 2 (l + r) 3 (l + 3r-3r 2 -r 3 



+ 4r 4 (l + r 2 

(l+r 2 )/2 



tan 



-i 



2r 



7T 



1-rV 2 
cfe 2t(x + r)xi(x — r/2) £(x) 2 



2r 4 (l-r) 2 lnr (35) 



— r 2 (l - r 2 )(l + 3r + 34r 2 + 3r 3 + r 4 ) 



+ r 4 (l +r) 



tan 



2r 
1 — r A 



71 



+ r 4 (2-r + 2r 2 )lnr (36) 
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l T2 



(l+r 2 )/2 



12 



dx 2t(x + r)x\{x + r/2) £(xy 
r 2 (l - r 2 )(l + 9r - 14r 2 + 9r 3 + r 4 



+ r 4 (3(1 + rf - 8r 



tan 



2r 
1 — r 2 



7T 



3r 5 lnr (37) 
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FIGURE CAPTIONS 
Fig. 1: Diagrams which contribute to M — > rrv^lu from (a) the SM W 
exchange, (b) charged Higgs exchange, (c) Wl-Wr mixing. 

Fig 2: j Ch+ as a function of r = m/M for £(x) given by the monopole 
approximation (|31|) with p = 1 (solid lines), p = 1.2 (dashed lines) and where 
£(x) 2 is naively divided out (dash-dot lines). The top, middle, and bottom 
sets of curves correspond to M -> mlu, M — > m* L li>, and M — > m*lv decays, 
respectively. 

Fig. 3: Pj L /C w + as a function of r. Notation is the same as in Fig. 2, 
with the top and bottom sets of curves corresponding to M — > m* T ^lv and 
M — > m* T ^lv decays, respectively. 
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